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ABSTRACT 


The  NASTRAN  Contour  Plotter,  a  group  of  subroutines  and  modifications  to 
the  NASTRAN  Plot  Module,  enable  contour  lines  to  be  superimposed  on  the  plot 
of  the  structural  model  or  on  an  outline  of  the  structural  model.  The  NASTRAN 
Contour  Plotter  can  be  incorporated  into  NASTRAN  Version  12.  Consistent  with 
the  NASTRAN  computer  program,  it  is  operational  on  the  IBM  360,  the  CDC  6000, 
and  the  UNIVAC  1108  computers  on  a  variety  of  plotter  hardware. 


ADMINISTRATIVE  INFORMATION 


The  work  was  performed  under  Subproject  SF  54  532  106, 

Task  15326,  sponsored  by  the  Naval  Ship  Systems  Command,  Code  0311. 


INTRODUCTION 

NASTRAN  (NAsa  STRuctural  ANalysis)program  is  a  general  purpose  digital 
computer  program  for  the  analysis  of  large  complex  structures.  NASTRAN  uses 
the  finite  element  structural  analysis  approach  which  involves  replacing  the 
actual,  continuous  structure  by  an  idealization  composed  of  small  elements. 

NASTRAN  currently  provides  a  plotting  capability  to  aid  in  tha  preparation 
of  the  structural  model  and  in  the  interpretation  of  the  computed  results. 

This  structure  plotting  capability  produces  both  deformed  and  undeformed 
geometric  projections  of  the  structural  model  with  or  without  displacement 
vectors-  at  the  nodal  points.  The  X-Y  plotting  capability  produces  graphs  of 
transient  response  and  frequency  response,  and  the  matrix  plotting  capability 
generates  topological  displays  of  matrices. 

The  contour  plotting  capability,  an  enhancement  to  the  structure  plotter, 
will  plot  lines  representing  equal  values  on  any  deformed  or  undeformed 
geometric  projection  of  the  structural  model.  It  was  developed  to  aid  in  the 
interpretation  of  various  stress  and  displacement  results. 

The  contour  plotting  capabilities  are  described  in  terms  of  the  data  card 
specifications  which  indicate  the  type  of  contour  plot  and  direct  the  contour 
plotter  to  superimpose  contour  lines  on  the  structural  model.  These  data  cards 
are  an  addition  to  the  structure  plot  request  packet  used  to  define  and 
control  the  execution  of  the  structural  model  plot  (ref.  1). 

Adding  the  contour  plotting  capability  to  the  structure  plot  module 
necessitated  certain  modifications  to  existing  subroutines  within  the  module 
and  to  the  plot  DMAP  (Direct  Matrix  Abstraction  Program)  statement  in  the  all 
rigid  formats.  These  modifications  and  additions  are  described  in  the  format 
of  The  NASTRAN  Programmer's  Manual  (ref.  2). 
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USING  THE  CONTOUR  PLOTTING  CAPABILITY 


The  Contour  Plot  Definition  Card 


The  contour  plot  definition  card  (fig.  1)  is  used  to  specify  the  type  of 
contour  plot  and  the  contour  values  to  be  plotted.  The  card  is  optional  since 
all  parameters  may  be  defaulted.  The  card  follows  the  format  of  a  free-field 
plot  request  card  (ref.  3). 


C0NT0UR 


stress 

displacement 


EVEN  10 
LIST  cl,c2 
LEVEN  n 


cn 


Z1 

Z2 

MAX 

MID 


—  * 


Figure  1. 


Nine  types  of  stress  contour  plots  are  available: 


stress 


MAJPRIN  \ 
MINPRIN  J 
MAXSHEAR 
J  XN0RMAL  ) 
S  YN0RMAL  > 
ZN0RMAL  ) 
XY SHEAR  ) 
XZSHEAR  J 
YZSHEAR  ) 


-  Major  and  Minor  Principal  Stress 

-  Maximum  Shear  Stress 

-  X ,Y ,Z  Normal  Stresses 


-  Shear-XY,  -XZ ,  -YZ  . 


Four  types  of  displacement  contour  plots  are  available: 


displacement 


XDISPLAC 

YDISPLAC 

ZDISPLAC 


the  X,Y,Z  components  of  the 
displacement  vector 


MAGNITUD  -  the  magnitude  of  the  displacement  vector 


Default  for  this  parameter  is  MAJPRIN,  Major  Principal  Stress. 


The  contour  values  to  be  plotted  may  be  specified  by  supplying  the 
parameter  EVEN  and  the  number,  n,  of  contour  values,  or  by  supplying  LIST  and 
a  list  of  real  number  contour  values.  If  EVEN  is  specified,  the  contour 
plotter  will  calculate  n  contour  values  at  (n-1)  equal  intervals  over  the 
range  of  values  specified  by  the  user.  The  first  contour  value  will  be  the 
minimum  and  the  n-th  contour  value  the  maximum  of  the  values  for  the  current 
set  of  elements.  The  number  of  contour  values,  n,  must  be  less  than  or 
equal  to  50.  Default  for  this  parameter  is  "EVEN  10",  i.e.,  contour  values 
at  10  even  (i.e.,  equal)  intervals. 

Since  stress  may  be  calculated  at  two  Fibre  Distances,  the  Fibre  Distance 
may  be  specified  by  designating  Z1  (Fibre  Distance  1)  or  Z2  (Fibre  Distance  2) . 
In  addition,  MAX,  the  maximum  of  Z1  and  Z2,  or  MID,  the  average  of  Z1  and  Z2 , 
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may  be  selected.  In  the  set  of  elements  for  which  a  contour  plot  can  be 
executed  (see  The  PLOT  Execution  Card),  the  average  of  Z1  and  Z2 ,  MID,  is 
applicable  only  to  the  CTR1A2  and  CQUAD2  elements.  The  default  is  Zl. 

The  Normal  Stresses  are  directional  and  are  calculated  in  the  elements' 
local  coordinate  system.  If  COMMON  DIRECTION  is  specified,  the  contour  plotter 
will  transform  the  Normal  Stresses  and  the  Shear  Stresses  to  a  common 
(currently  the  basic)  coordinate  system.  If  LOCAL  DIRECTION  is  specified, 
the  contour  plotter  will  leave  these  stresses  in  the  global  coordinate  system. 
Note  that  the  Normal  Z  Stress,  and  the  Shear-XZ  and  Shear-YZ  are  assumed  to  be 
zero  in  the  element's  local  coordinate  system.  Default  for  this  parameter  is 
the  COMMON  DIRECTION. 


The  PLOT  Execution  Card 


The  PLOT  Execution  Card  (figure  2)  has  been  altered  to  include  the 
CONTOUR  and  OUTLINE  execution  commands  (ref.  4). 


{STATIC  )  j 

MODAL  >  DEFORMATION  il, 12, THRU  i3,i4,  etc.  > 
TRANSIENT)  J 

[■/RANGE  XI,  A 2)1 
JjTIME  tl,  1 2/ J 

[CONTOUR  il ,i2  THRU  13, i4,  etc.]  >  [{^TtlJ ’t^  }]  > 

K  SYMMETRY  l  ~\ 

antisymmetry}  wJ 
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(DENSITY 


|  pj  >  ^SYMBOLS  m[,n]J>|i 


(GRID  POINTS 
ELEMENTS 
BOTH 
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OUTLINE 
VECTOR  v 
SHAPE,  VECTOR  v 


Figure  2. 


The  CONTOUR  command  will  superimpose  lines  representing  equal  value  on 
the  structural  model  or  on  the  outline  of  the  structural  model.  The  nonzero 
integers  following  the  CONTOUR  command  refer  to  subcases  to  be  plotted. 
Default  is  the  subcase  following  the  DEFORMATION  parameter  or,  if  that 
parameter  is  not  specified,  all  subcases.  In  addition,  the  RANGE  or  TIME  may 
be  specified.  If  both  DEFORMATION  and  CONTOUR  are  specified,  the  contour 
plotter  will  take  the  subcase,  range,  or  time  parameters  following 
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DEFORMATION  as  applicable  to  the  contour  plot  as  well  as  to  the  deformed  plot. 

If  a  stress  contour  plot  is  requested,  contour  lines  will  be  plotted  on 
the  following  elements  in  SET  j  appropriate  to  the  type  of  stress  contour 
plot  requested  (ref.  5):  CSHEAR,  CTR1A1 ,  CTR1A2 ,  CQUAD1 ,  CQUAD2 ,  CTRMEM , 
CQDMEM. 

If  a  displacement  contour  plot  is  requested,  contour  lines  will  be  plotted 
on  all  two-dimensional  elements  plotted  by  the  structure  plotter.  The  type  of 
deformation  must  be  specified  for  a  displacement  contour  plot. 

The  contour  lines  are  labeled  with  integers  indicating  the  contour  value. 
The  integers  are  listed  with  their  associated  contour  values  in  the  form  of  a 
key  in  the  Plot  Module  Message  section  of  the  printed  output. 

The  OUTLINE  command  has  been  added  as  an  alternative  to  the  SHAPE 
command.  OUTLINE  will  plot  connecting  lines  between  those  grid  points  which 
lie  on  the  boundary  of  the  structural  model. 

If  DEFORMATION  is  specified  in  an  addition  to  OUTLINE,  the  OUTLINE  will 
reflect  the  deformed  structure.  If  a  "0"  is  added  to  the  deformation  subcase 
string,  the  undeformed  outline  will  be  plotted. 

The  OUTLINE  command  may  be  specified  only  in  addition  to  the  CONTOUR 
command . 


CONTOUR  PLOT  PROGRAM  FUNDAMENTALS 
Modifications  and  Additions  to  the  Functional  Module  Plot 

The  addition  of  the  contour  plotter  enlarges  the  capabilities  of  the 
Functional  Module  Plot  (ref.  6)  to  include  the  capability  of  superimposing 
contour  lines  on  structural  shapes  using  a  variety  of  different  plotters.  In 
order  to  incorporate  the  contour  plotter,  two  data  blocks  were  added  to  the 
input  string  of  data  blocks  brought  into  the  Plot  Module  and  modifications 
were  made  to  the  existing  subroutines  in  order  to  add  the  new  subroutines  to 
the  Plot  Module.  A  page  of  printed  output  was  also  devised  to  give  information 
on  the  contour  plot. 

The  data  blocks  ECPT  and  0ES1  were  added  as  additional  input  blocks  to 
the  DMAP  calling  sequence.  The  DMAP  calling  sequence  (ref.  7)  becomes: 

PL0T  PLTPAR,  GPSETS ,  ELSETS  ,  CASECC,  BGPDT ,  EQEXIN ,  SIL,  ECPT, 

PLTDSP1 ,  PLTDSP2 ,  0ES1/PL0TX/V ,N ,NGP/V ,N ,LSIL/V ,N ,PLTFLG/ 

V  ,N ,PLTNUM  $ 

ECPT  is  the  Element  Connection  and  Properties  Table  and  0ES1  is  the  Output 
Element  Stress  Table.  If  the  ECPT  data  block  is  purged,  contour  plots  will 
not  be  drawn.  If  the  0ES1  da'ta  block  is  purged,  the  stress  contour  plots 
will  not  be  drawn. 
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Modifications  were  made  to  existing  plot  module  subroutines  so  that  the 
contour  plot  definition  card  could  be  detected  and  so  that  the  CONTOUR  and 
OUTLINE  commands  could  be  executed.  When  the  plot  module  is  ready  to  generate 
the  plot  (subroutine  DRAW),  all  the  necessary  parameters  have  been  set  and 
any  necessary  rotations  have  been  made  to  generate  the. contour  plot. 

The  subroutine  CONTOR  is  called  by  the  subroutine  DRAW  to  plot  contour 
lines.  If  a  stress  contour  plot  was  requested,  CONTOR  requests  that  the  plot 
set  be  ordered  and  that  contour  values  be  found  for  each  element.  The 
contour  lines  are  then  located  on  the  structural  model,  drawn  and  labeled. 

If  a  displacement  contour  plot  was  requested,  the  contour  lines  are  drawn  by 
the  subroutine  DISPLA  which  returns  to  the  subroutine  CONTOR  to  plot  the 
labels . 

If  an  OUTLINE  of  the  structural  model  was  requested,  the  subroutine  DRAW 
calls  the  subroutine  BORDER  to  generate  this  part  of  the  plot. 

The  subroutine  DRAW,  then,  (ref.  8),  has  been  altered  so  that  subroutines 
CONTOR  and  BORDER  can  be  called  to  generate  contour  lines  and  the  structural 
outline.  The  subroutine  PLOT  (ref.  9)  has  also  been  altered  so  that  a  call 
to  ROTATE  can  be  made  to  transform  certain  stresses  to  a  common  coordinate 
system. 

The  seven  subroutines  comprising  the  contour  plotting  capability  are 
described  below. 

(1)  Subroutine  Name:  R0TATE 

Purpose:  to  transform  the  Normal  Stresses  and  the  Shear  Stresses  to 

a  common  coordinate  system. 

Calling  Sequence:  CALL  R0TATE(ELSETS ,  GPLST,  X) 
where  ELSETS  is  the  data  block  containing  element  connection  data, 

GPLST  is  a  list  of  indices  defining  the  set  of  grid  points  associated 
with  the  elements,  and 

X  contains  the  corresponding  coordinates  of  the  grid  points. 

Common  Blocks : 

/XXPARM/  (ref.  10)  (see  also  Additions  to  NASTRAN  Tables, 

page  10  ) . 

Method:  The  subroutine  transforms  the  Normal  Stresses  and  Shear  Stresses 
from  the  local  coordinate  system  to  a  common  coordinate  system  using  the 
following  expression: 


where  [T] 


[T]  =  [6]  [T']  te] 


xy 


xz 


yx 


zx 


y  yz 

T  O 

zy  z 


is  a  symmetric  matrix. 
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are  the  Normal  X,Y,  and  Z  Stresses  respectively, 


t  ,  t  ,  t  are  the  Shear-XY,  -XZ,  and  -YZ  Stresses 
xy  xz  yz  ,  , 

respectively,  and 

[0]  is  the  rotation  matrix  which  has  as  its  components  the 
direction  cosines  between  the  local  and  common 
coordinate  systems. 

Since  all  elements  considered  by  the  contour  plotter  are  plate  elements, 

the  elements  are  all  assumed  to  lie  in  the  local  X,Y  plane  with  the  X-axis 

passing  from  grid  point  one  to  grid  point  two  (ref.  11).  The  subroutine 

makes  use  of  the  assumption  that  o=t  =  t  =0in  the  local  coordinate 

z  xz  yz 

system.  J 

A  new  local  element  stress  table  (NEWOES)  is  compiled  by  the  subroutine 
ROTATE  to  communicate  the  transformed  stresses  to  the  subroutine  CONTOR. 


Additional  Subroutines  Required:  None. 


(2)  Subroutine  Name:  C0NT0R 


Purpose:  to  plot  and  label  contour  lines. 


Calling  Sequence:  CALL  C0NT0R(GPLST ,X ,S ,U ,PEN .DEF0RM) 
where  GPLST  is  a  list  of  indices  defining  the  set  of  grid  points, 

X  contains  the  corresponding  coordinates  of  the  undeformed  points, 

S  contains  the  corresponding  coordinates  of  the  deformed  points, 

U  contains  the  displacement  coordinates, 

PEN  is  the  line  density  or  pen  number,  and 


DEFORM  = 


0  if  an  undeformed  structural  model  was  requested 
1  if  a  deformed  structural  model  was  requested 


Common  Blocks : 

/XXPARM/  (ref. 12)  (see  also  Additions  to  NASTRAN  Tables, 

page  10  ) . 


/PLTDAT/  (ref. 13) 


/DRWDAT/  (ref. 12)  (see  also  Additions  to  NASTRAN  Tables, 

page  10  ) . 


Method:  If  the  contour  plot  is  to  be  a  displacement  contour  plot, 

CONTOR  calls  the  subroutine  DISPLA  to  plot  the  contour  lines.  If  a  stress 
contour  plot  was  requested,  CONTOR  calls  the  subroutine  CREATE  to  prepare  a 
table  containing  the  contour  value  and  the  centroid  for  each  element  in  the 
plot  set.  The  subroutine  ORDER  is  then  called,  if  necessary,  to  provide  a 
table  of  all  grid  points  internal  to  the  structural  shape  and  all  elements 
that  share  a  common  grid  point  in  the  current  plot  set. 


Once  these  tables  are  available,  CONTOR  estimates  a  contour  value 
for  each  internal  grid  point  using  a  weighted  average  of  distance  from  the 
centroids  and  contour  values  of  surrounding  elements .  Contour  lines  are 
found  within  a  triangle  formed  by  connecting  an  internal  grid  point  and  two 
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adjacent  elements'  centroids.  Contour  lines  are  plotted  with  the  assumption 
that  the  contour  values  vary  linearly  between  vertices  of  this  triangle. 


Labels  are  placed  on  the  completed  contour  plot  and  the  subroutine 
PLTOPR  is  called  to  print  the  contour  plot  key  as  part  of  the  plot  messages. 

Additional  Subroutines  Required:  ORDER,  CREATE,  DISPLA. 


(3)  Subroutine  Name:  CREATE 

Purpose:  to  extract  the  stress  contour  values  from  the  appropriate  data 

block  and  calculate  the  centroid  of  each  element  in  the  plot  set. 


Calling  Sequence:  CALL  CREATE (GPLST ,X ,U , DEF0RM , C0NMIN , C0NMAX , ELMTID , ST0RE) 
where  GPLST  is  a  list  of  indices  defining  the  set  of  grid  points, 

X,U  contain  the  corresponding  undeformed  and  deformed  coordinates 
respectively,  of  the  grid  points, 


DEFORM  = 


{! 


if  an  undeformed  structural  shape  was  requested, 
if  a  deformed  structural  shape  was  requested 


CONMIN.CONMAX  are  the  minimum  and  maximum  of  the  contour  values,  and 
ELMTID, STORE  are  temporary  storage  areas  for  the  element 

identification  number  and  the  contour  value  found  on 
the  stress  data  block. 


Common  Block: 

/XXPARM/  (ref. 14)  (see  also  Additions  to  NASTRAN  Tables,  page  10  ). 

Method:  After  determining  which  stress  data  block  contains  the  contour 
values  for  this  contour  plot,  the  appropriate  data  block  is  opened  and  the 
temporary  storage  area  is  filled.  The  contour  values  are  selected  from  the 
information  on  Fibre  Distance,  direction,  subcase,  and  if  applicable,  range 
or  time  specified  by  the  user.  Then  the  ELSETS  data  block  is  read  to  supply 
information  on  element  types,  elements  identification  numbers,  and  internal 
grid  point  numbers  of  elements  contained  in  the  current  set. 


The  subroutine  makes  two  assumptions  concerning  the  stress  data  block 
and  the  ELSETS  data  block.  It  assumes  that  element  types  are  ordered  in  the 
same  sequence  and  that  for  each  element  type,  element  information  is  ordered 
by  increasing  element  identification  number.  Element  types  for  which  contour 
lines  cannot  be  drawn  are  ignored  as  are  elements  not  included  in  the  current 
plot  set. 


Finally,  the  centroid  is  calculated  for  each  element  in  the  plot  set 
and  a  temporary  table  is  made  containing  the  element  identification  number, 
the  contour  value,  and  the  centroid. 

Additional  Subroutines  Required:  CENTRE  . 
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(4)  Subroutine  Name:  0RDER 

Purpose:  to  find  for  each  grid  point  in  the  plot  set,  all  elements 

'which  contain  that  grid  point  and  to  order  those  elements  surrounding  interior 
grid  points. 

Calling  Sequence:  CALL  0RDER(GPLST ,  ID,  REST) 
where  GPLST  is  a  list  of  indices  defining  the  set  of  grid  points  and 
ID, REST  are  temporary  storage  areas  for  element  and  grid  point 
information . 

Common  Blocks:  None. 

Method:  The  subroutine  ORDER  extracts  from  the  ECPT  (Element 
Connections  and  Properties  Table)  data  block  the  list  of  elements  connected 
to  each  grid  point  for  each  grid  point  in  the  plot  set.  Since  the  grid 
points  numbers  given  in  the  ECPT  data  block  are  in  SIL  (Scalar  Index  List) 
form,  ORDER  uses  the  SIL  table  to  find  the  internal  grid  point  numbers 
associated  with  each  element.  For  each  grid  point  in  the  plot  set,  ORDER 
finds  all  two  dimensional  elements  in  the  plot  set  containing  that  (pivot) 
grid  point.  The  elements  are  then  ordered  so  that  adjacent  elements  contain 
a  second  common  grid  point.  A  flag  is  set  indicating  whether  the  pivot  grid 
point  is  interior  to  the  structural  shape  or  on  the  boundary  of  the  structural 
shape. 


A  temporary  table  is  compiled  containing  one  record  for  each  grid 
point  in  the  current  plot  set.  Each  record  contains  a  flag  indicating 
whether  the  grid  point  is  interior  to  the  structural  model  or  on  the  boundary 
of  the  structural  model,  and  a  list  of  elements  connected  to  the  grid  point. 

For  boundary  grid  points,  the  table  contains  two  grid  points  for  each  element 
that  are  connected  by  element  boundary  lines  to  that  grid  point. 

Additional  Subroutines  Required:  None. 

(5)  Subroutine  Name:  CENTRE 

Purpose:  to  find  the  intersection  of  two  lines. 

Calling  Sequence:  CALL  CENTRE (XI ,Y1 ,X2 ,Y2 ,X3 ,Y3 ,X4 ,Y4 , CENTER) , 

RETURNS (RETURN1) 

where  Xl,Yl  and  X3,Y3  are  the  end  points  of  one  line, 

X2,Y2  and  X4,Y4  are  the  end  points  of  the  second  line,  and 
CENTER  is  returned  as  the  point  of  intersection. 

Common  Blocks:  None. 

Method:  The  subroutine  CENTRE  is  used  in  the  calculation  of  the  centroid 
of  an  element.  The  nonstandard  return  is  taken  in  the  event  that  either  of  the 
lines  has  infinite  slope.  It  is  assumed  that  the  two  lines  are  not  parallel. 

Additional  Subroutines  Required:  None. 
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(6)  Subroutine  Name:  DISPLA 

Purpose:  to  draw  displacement  contour  lines. 


Calling  Sequence:  CALL  DISPLA (GPLST ,X , S , U , PEN , DEFORM , LABEL , PT) 
where  GPLST  is  a  list  of  indices  defining  the  set  of  grid  points, 
X,S  contain  the  corresponding  undeformed  and  deformed 
coordinates  respectively,  of  the  grid  points. 

U  contains  the  grid  point  displacement  coordinates, 

PEN  is  the  line  density  or  pen  number, 


DEFORM  = 


0 

1 


if  an  undeformed  structural  shape  was  requested 
if  a  deformed  structural  shape  was  requested  , 


LABEL  contains  the  contour  line  labeling  information,  and 
PT  contains  the  coordinates  of  the  triangle  formed  to  calculate 

contour  lines. 


Common  Blocks : 

/XXPARM/  (ref. 15)  (see  also  Additions  to  NASTRAN  Tables,  page  10  ). 

/PLTDAT/  (ref. 15)  (see  also  Additions  to  NASTRAN  Tables,  page  10  ). 

Method:  The  subroutine  DISPLA  uses  the  ELSETS  data  block  to  extract 
the  two  dimensional  elements  in  the  plot  set.  Each  element  is  then  broken 
down  to  one  or  more  triangles  and  contour  lines  are  drawn  within  the 
triangles  with  the  assumption  that  the  contour  values  vary  linearly  between 
grid  points. 


It  is  assumed  in  the  DISPLA  subroutine  that  the  user's  coordinate 
system  is  a  right-handed  coordinate  system  with  coordinate  axes  X,Y,Z.  If 
the  AXES  card  (ref.  16)  is  used  to  change  the  orientation  of  the  structural 
model,  it  will  have  no  effect  on  the  component  of  displacement  used  to  draw 
the  contour  plot. 


Labeling  information  is  maintained  on  a  temporary  file  to  be  added 
to  the  plot  by  the  subroutine  CONTOR. 

Additional  Subroutines  Requested:  None. 


(7)  Subroutine  Name:  B0RDER 


Purpose:  to  draw  the  structural  outline. 


Calling  Sequence:  CALL  B0RDER(GPLST,  X,  U,  DEF0RM) 
where  GPLST  is  a  list  of  indices  defining  the  set  of  grid  points, 

X,U  contain  the  corresponding  undeformed  and  deformed  coordinates, 
respectively,  of  the  grid  points,  and 


DEFORM 


-{! 


if  an  undeformed  structural  shape  was  requested 
if  a  deformed  structural  shape  was  requested. 


Common  Blocks  :  None . 


Method:  The  subroutine  B0RDER  uses  the  table  created  by  the  subroutine 
0RDER.  BORDER  extracts  from  the  table  those  grid  points  that  lie  on  the  plot 
sets  boundary.  All  grid  points  connected  to  this  grid  point  by  element 
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boundary  lines  are  sorted  to  find  the  two  element  boundary  lines  which  form 
the  structural  boundary.  The  BORDER  subroutine  then  draws  two  connecting 
lines  for  each  grid  point  so  that  the  structural  outline  will  be  of  density 
two . 


Additional  Subroutines  Required:  None. 


Additions  to  the  NASTRAN  Tables 


The  NASTRAN  Functional  Module  Plot  uses  labeled  common  to  communicate 
information  concerning  plot  requests,  plot  calculation  parameters,  and  plotter 
hardware  characteristics  among  subroutines.  Additional  parameters  have  been 
added  to  two  of  the  labeled  common  blocks  to  pass  contour  plotter  information. 
These  additional  parameters  are  described  here  (ref.  17). 

(1)  C0MM0N / DRWDAT / SET , . . . , VECT0R , C0N , EDGE 

where  C0N  is  an  integer  indicating  the  contour  plot  option,  and 
EDGE  is  an  integer  indicating  the  outline  option. 


(2)  C0MM0N/XXPARM/PBUFSZ , . . . ,XY(11,3) ,NCNTR,CNTR(50) , ICNTVL , WHERE , DIRECT , 
SUBCAS , FLAG .VALUE , SET 


In  the  following  descriptions,  the  value  in  parenthesis  to  the  right  of 
the  variable  name  is  the  default  value,  and  the  letter  in  parenthesis  to  the 
right  of  the  explanation  pertains  to  the  type  of  variable  (I  is  integer,  and 
R  is  real) . 


Contour  Plotting  Data 


NCNTR(IO) 

CNTR 

ICNTVL (3) 
WHERE (1) 
DIRECT (2) 
SUBCAS 
FLAG (0.0) 
VALUE 
SET 


Number  of  contour  values  to  be  plotted  (I) 
Contour  values  (R) 

Type  of  contour  plot  (I) 

Fibre  Distance  of  stress  contour  value  (I) 
Stress  vector  direction  option  (I) 

Current  subcase  (I) 

Data  identification  indicator  (R) 

Current  eigenvalue  or  time  step  (R) 

Last  plot  set  processed  by  contour  plotter  (I) 


Local  Table  Descriptions 


Three  tables  are  formed  for  the  use  of  the  contour  plot  subroutines.  The 
tables  are  compiled  and  used  within  the  structure  plot  module.  The 
description  of  each  table  follows  the  file  name. 

1.  SCR1,  scratch  file  1  -  compiled  by  the  subroutine  CREATE  to  be  used  by 
the  subroutine  CONTOR.  The  header  record  is  written  by  the  subroutine 
GOPEN.  The  table  is  compiled  from  information  contained  in  the  ELSETS  data 
block  (ref.  18)  and  the  0ES1  data  block  (ref.  19)  or  the  NEWOES  table  (see 
below) . 
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Table  Format: 


Record  Word  Type 

0 

111 

2  I 

3  R 

4,5  R 

(4m+2)  I 


Item 

Header  Record 


Element  Type 


Element 

Identification 
Numb  er 

Stress  Value 
Element  Centroid 


repeated 
^  for  each 
element 
in  plot 


repeat  for 
each  element 
type  in 
contour  type 
set 


2 


End-of-File 


Note:  m  =  the  number  of  elements  of  current  type  in  the  plot  set. 


2.  SCR2,  scratch  file  2  -  compiled  by  the  subroutine  ORDER  to  be  used  by 
the  subroutines  CONTOR  and  BORDER.  The  header  record  is  written  by  the 
subroutine  GOPEN.  The  table  is  compiled  from  information  contained  in  the 
ECPT  data  block  (ref.  20)  and  the  SIL  table  (ref.  21). 


Table  Format: 


Record  Word 
0 

1  1 
2 
3 


4 


Type 

I 

I 

I 

I 


Item 

Header  Record 


Flag 

Number  of  elements,  m. 


Grid  point  internal 
number  ( 

repeated  for  each 
element  containing 
the  grid  point. 

Contents  dependent 
on  Flag.  _ / 


repeated  for 
each  grid  point 
in  plot  set 


Notes : 

1)  n  =  the  number  of  grid  points  in  the  plot  set. 

2)  m  =  the  number  of  elements  in  current  record. 

3)  Flag  =  -1: 

Word  Type  Item 

4  I  Element  identification  number 


(3+m)  I 
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4) 


Flag  =  -2: 
Word 


4 

5 

6 


Type  Item 

I  Element  identification  number 

I  Grid  point  internal  number 

I  Grid  point  internal  number 


(3m+l)  I 
(3m+2)  I 
(3m+3 )  I 


Element  identification  number 
Grid  point  internal  number 
Grid  point  internal  number 


3.  NEWOES  -  compiled  by  the  subroutine  ROTATE  to  be  used  by  the  subroutine 
CREATE.  There  is  no  header  record.  The  table  is  compiled  using  stress 
values  contained  in  the  0ES1  data  block  (ref.  22). 


The  first  record  in  each  set  of  two  records  on  NEWOES  is  a  copy  of 
the  corresponding  record  of  the  0ES1  data  block.  The  second  record  of  each 
set  contains  changes  in  content  for  the  following  element  types: 

Element  types  CTRIA1,  CQUAD1,  CTRIA2 ,  CQUAD2 

Word  Contents 

2  Z1 

3  Normal  X  Stress  at  Z1 

4  Normal  Y  Stress  at  Z1 

5  Normal  Z  Stress  at  Z1 

6  Shear  XY  at  Z1 

7  Shear  XZ  at  Z1 

8  Shear  YZ  at  Z1 

9  Maximum  Shear  at  Z1 

10  Z2 

11  Normal  X  Stress  at  Z2 

12  Normal  Y  Stress  at  Z2 

13  Normal  Z  Stress  at  Z2 

14  Shear  XY  at  Z2 

15  Shear  XZ  at  Z2 

16  Shear  YZ  at  Z2 

17  Maximum  Shear  at  Z2  ' 

Element  types  CTRMEM,  CQDMEM 


Word  Contents 


2 

3 

4 

5 

6 

7 

8 


Normal  X  Stress 
Normal  Y  Stress 
Normal  Z . Stress 
Shear  XY 
Shear  XZ 
Shear  YZ 
Maximum  Shear 
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APPENDIX  -  SAMPLE  OUTPUT 


A  square  plate  with  built-in  edges  and  concentrated  transverse 
loads  was  used  to  generate  sample  contour  plots.  There  are  100 
elements  in  the  structural  model. 

For  each  output  example,  the  contour  plot  (see  Figures  3,4,5) 
is  shown  below  a  sketch  of  the  structural  model.  The  loads  are 
indicated  on  the  sketch  by  X's  on  the  appropriate  grid  points.  The 
integers  used  to  label  contour  lines  refer  to  a  key  of  contour  values 
in  the  Contour  Plot  Message  Section  of  the  printed  output  (see 
Figure  6) . 

A  list  of  the  Structural  Plot  Data  Cards  used  to  generate  the 
sample  contour  plots  follows: 


0UTPUT(PL0T) 

SET  1  INCLUDE  ALL 
MAXIMUM  DEF0RMATI0N  1.0 
FIND 

PL0T  C0NT0UR  2  (see  Figure  3) 

PL0T  CONTOUR  2  0UTLINE  (see  Figure  3) 

C0NT0UR  MINPRIN 

PL0T  STATIC  DEF0RMATI0N  C0NT0UR  1  0UTLINE  (see  Figure  4) 
C0NT0UR  MAXSHEAR 

PL0T  C0NT0UR  2  0UTLINE  (see  Figure  5) 
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Figure  6.  Sample  Output  -  Contour  Plot  Message  Page 
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